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Fermion dark matter (DM) interacting with the standard model through a Higgs portal requires 
non-renormalizable operators, signaling the presence of new mediator states at the electroweak scale. 
Collider signatures that involve the mediators are a powerful tool to experimentally probe the Higgs 
portal interactions, providing complementary information to strong constraints set by direct DM 
detection searches. Indirect detection experiments are less sensitive to this scenario. We investigate 
the collider reach for the mediators using three minimal renormalizable models as examples, and 
requiring the fermion DM to be a thermal relic. The Large Hadron Collider in its high-energy, high- 
luminosity phase can probe most scenarios if DM is lighter than about 200 GeV. Beyond this scale, 
future high-energy experiments such as an electron-positron collider or a 100-TeV proton-proton 
collider, combined with future direct detection experiments, are indispensable to conclusively test 
these models. 


I. INTRODUCTION 

Dark matter (DM) can couple to a gauge invariant operator W H, where H is the standard model (SM) 
Higgs field. Since WH is the lowest-dimensional Lorentz- and gauge-invariant operator in the SM, such Higgs 
portal couplings to DM could dominate the interactions between the visible sector and the dark sector [iHin!. 
If DM is a scalar, (p, the Higgs portal operator, {(p^ip){H^H), is renormalizable and the theory is in principle 
UV-complete. In contrast, if DM is a fermion, y, the Higgs portal interactions are of mass dimension five, 

^EFT = + '^{xXbX){H^ H). (1) 

This necessarily implies the existence of new states at an energy scale A, the mediators between the Higgs field 
and the DM. Integrating out the mediators gives the above interactions in an Effective Field Theory (EFT). 
The presence of mediators offers new opportunities to experimentally search for the Higgs portal: Instead of 
searching directly for DM, one can search for the mediators instead. 

The mediators in general fall into two categories: i) mediators that are electroweak singlets, or ii) mediators 
that are charged under the electroweak gauge group. In this paper we explore the phenomenology of both 
categories using three minimal renormalizable models of the Higgs portal (henceforth called UV completions). 
As a representative of the hrst category we choose a model where both the DM and the mediator are electroweak 
singlets {the singlet-singlet model). For the second category we consider two examples, a model where DM 
can be an electroweak singlet {the singlet-doublet model), and a model where DM is necessarily part of an 
electroweak multiplet {the doublet-triplet model). 

For each model we compute the predicted relic density from thermal freeze-out and determine the parameter 
space that is consistent with the observed dark matter relic density. The viable parameter regions are con¬ 
fronted with bounds from direct detection, indirect detection and collider experiments. Interestingly, we find 
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that the combination of these constraints often requires the mediators to have masses of similar magnitude as 
the DM fermion, This means that the EFT Lagrangian in Eq. 0 is in general not a good description for 
the computation of the relic density and the collider phenomenology. For these observables, we therefore treat 
all fields in the dark sector as dynamical degrees of freedom, i.e. we “integrate in” the Higgs portal. 

Finally, we analyze how the UV Higgs portal completions can be probed by future direct detection and 
collider experiments. In particular, we consider the upcoming 14-TeV run of the Large Hadron Collider (LHC) 
with up to 3000 fb“^ luminosity, a future e“'"e“ collider with a center-of-mass energy of up to ^/s = 1 TeV, 
and a future pp collider with ^/s = 100 TeV. 

The paper is organized as follows. In Sec. [H] we introduce the three minimal UV completions of the fermion 
DM Higgs portal, in Sec. m we calculate the respective thermal relic densities, and in Sec. |IV| we estimate 
the DM direct detection rates. Sec. |V] deals with the constraints following from Higgs decays. In Sec. VI 


we 


combine all the above constraints with the expected sensitivity of the LHC 14-TeV run and future colliders 
on the mediators. Appendix contains the analytic formulae for loop-induced DM-Higgs interactions in a 
version of the doublet-triplet model that give the dominant direct detection signal. 


II. THE MODELS 

We start by introducing three minimal UV completions of the Higgs portal models with fermionic DM: 
(A) the singlet-singlet model, (B) the singlet-doublet model, and (C) the doublet-triplet model. All three 
examples are treated as simplified models rather than full theories, i.e. we do not consider issues such as 
anomaly cancellation and stability under renormalization group running to high scales. 


A. The Singlet-Singlet Model 

The dark sector is assumed to consist of DM, which is a Z 2 -odd SM-singlet Majorana fermion, x, and of a 
mediator, which is a Z 2 -even real singlet scalar, fs- All SM fields are Z 2 even. This model has already been 
discussed at length in the literature, see for example Refs. [iiHn]. In this paper we update the limits on the 
model and show projections for the 14-TeV LHC and a future 100-TeV pp collider. 

The relevant terms for interactions between the dark sector and the SM in the Lagrangian are 

£s D - - Ao|iL|" - - Uo(<(>s) - - A'o|LIp4 

1 1 ■ 

- 2^x,oXX - ^yoXX^ks - |2/5.oX75X<(’s, 

where H is the SM Higgs field and Vo{(f)s) contains cubic and quartic (j)s self-interactions. Here we use the four- 
component notation for the Majorana fermion y. The scalar fields H and 4>s acquire the vacuum expectation 
values (vevs) v = 246 GeV and vs, respectively. The vev v breaks the electroweak symmetry, while vs does 
not. We can therefore write 


H = {G+,{v + h + iG°)/V2y, cl)s = vs + S, (3) 

with the would-be Goldstone bosons eaten by the longitudinal components of the IU+ and Z bosons, 

respectively. In terms of h and S the interaction Lagrangian is 

£s D ^2+ ^ + - V{S) - (p's + X'S^){vh + 

- h^xxx - \yxxS - ^VbXibxS, 

where -I- yovs is the DM mass, y = yo and 2/5 = 2 / 5,0 are the parity-conserving and parity-violating 

Yukawa couplings, respectively, and m| = mg q + Xqv'^ + Vq {vs) is the singlet mass squared. The Higgs-singlet 

mixing parameter is given by p' = + 2 AqVs, while = Pq + PqVs + X'^Vg, and A^ ^ = Aq^ We also define 

V{S) = psS^ + XgS'^ that contains triple and quartic singlet scalar interactions. 
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The DM state x interacts with the SM through Higgs-singlet mixing. For /r' 7^ 0, the mass eigenstates hi ^2 
are admixtures of h and S, 


hi — Cah SaSj 


/12 — Sah CaS. 


(5) 


Here we use the abbreviations Ca = cos a, Sq = sin a. The mixing angle a and the masses of hi ^2 are given by 


2^1/ m|+^\ 

“ 2 V /’ 



2^^ =F , 


( 6 ) 


with the mass splitting (Am^)^ = (m| + 2/i^)^ + The state hi is the observed SM-like scalar with 

mass mti = 125.09 ± 0.24 GeV [H], whereas ft.2 is mostly singlet-like. The mixing angle a is constrained by 
the measured Higgs production and decay rates and, depending on the mass mh 2 , by the non-observation of 
a second Higgs scalar at the LHC, as we show in Sec. |VI A[ In particular, Higgs-singlet mixing leads to the 
DM-SM interactions (here and below C denotes Lagrangian in the mass eigenstate basis) 


^ A -\y{ 


, xxhl + Ca xxh^) - §2/5 {-Sa XlbXhl + Ca X75X^2) ffhl + Sa ///12) , (7) 


where y/ is the SM Yukawa coupling of the fermion /. The scalar interactions also govern direct detection 
signatures of spin-independent DM scattering off nuclei. Since the relevant effective coupling XXQQ is ~ sin(2a) 
(see Tab.|T|, direct detection experiments set strong constraints on the mixing angle a. (see Figs. |^and|^. 


B. The Singlet-Doublet Model 

In this model the dark sector consists of two fermion fields, xd and xs-, transforming under the SM elec- 
troweak group SU{2)l x U{1)y as 


Xx5~(2,1/2), xs~(1,0). (8) 

The field xd = (Xd^Xd) i® a doublet of Dirac fermions with vector-like gauge interactions, while the singlet 
Xs can be either a Dirac or a Majorana fermion. We discuss both of these possibilities. In the dark sector 
we impose a Z2 symmetry, xd,s ~Xd,Si under which the SM fermions are even. This forbids mixing with 
neutrinos and makes the lighter of the two mass eigenstates in the dark sector stable. 

a. Dirac singlet fermion. We first discuss the case where xs is a Dirac fermion. The particle content of 
this model consists of two neutral Dirac fermions, xs and x^: and a charged Dirac fermion x~d- The relevant 
terms in the Lagrangian are 


Cm A -moXDXD - msxsxs - {yXoXsH -k h.c.). (9) 

Without loss of generality, the Yukawa coupling y can be chosen real by redefining the complex phase of the 
Xd- After electroweak symmetry breaking (EWSB) this term introduces the mixing between xs and x%- For 
later convenience, we define the mixing angle 9a generally as 

sin^ ffa = ^ ^1 -k ^ ^ ’^Th = {ms - moY + a{yvY . (10) 

Here a = 2, and the heavy and light mass eigenstates, Xh and x?j are given by 

X° = COS02XS+ sin6»2x2), X? = - sin 6»2 xs + cos 6*2 Xd , (11) 

with the corresponding mass eigenvalues 

+ TOs ± Am2). 


( 12 ) 
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The charged state Xd ^ mass m'^ = mo- In the mass eigenstate basis, the interactions of the neutral 
fermions with the Z and Higgs bosons read 


£ D - 


2c\y - 


cos^ 6»2 X?7^x? + sin^ 62 xll^xl + ^ sin(202) (x°7^X? + X?7^X°) 


(13) 


y 


sin(202) {thXl 


X°X°) + cos(26l2) (xftX° + X°X°) 


h. 


It is interesting to consider two parameter limits of this model. For |m£) — ms\ «C \yv\^ xs and x% 
maximally mixed, 62 ~ 7r/4. The two neutral mass eigenstates are split by — mf « v^lynl, while the 
charged state has a mass = ttid ~ + m°)l2, and thus lies in between the two neutral states. Because 

of the large mixing, the DM coupling to the Z boson is unsuppressed, so that direct DM detection searches 
exclude this possibility. 

If |tod — ms\ S> yv, Xh i® significantly heavier than x?- For mu > ms, the DM state is mostly a singlet, 
X? - -XS + O'xd - \yv\/V^{mD — ms). The coupling of DM to the Z boson is thus suppressed, and 

the model is allowed by direct DM searches. In contrast, if ms > mo, DM is mostly a doublet, x? Xd^ with 
unsuppressed couplings to the Z boson. This possibility is therefore excluded by direct DM searches. 


b. Majorana singlet fermion. The second possibility is that xs i® & Majorana fermion. This scenario 
corresponds to the bino-higgsino system (with decoupled wino) in the Minimal Supersymmetric Standard 
Model (MSSM) for tan/3 = 1 and y = g'l\/2. Some phenomenology of the singlet-doublet Majorana model 
has also been explored in Refs. [T9H23]^. We will slightly abuse the notation, such that in this subsection xs 
denotes a two-component Weyl fermion, while xv xfh Iwo left-handed Weyl fermions forming a Dirac 
fermion. They transform under SU{2)l x U{1)y as xd ~ (2,1/2) and Xri ~ (2,—1/2). The translation to 
four-component notation is given by (xd , cXd ) Xd , where the final xd is the Dirac fermion in Eq. ([^ (here 
eij is the antisymmetric tensor in the SU{2)s space). In two-component notation, the relevant terms in the 
Lagrangian read 

£m £ moxh^XD - ^msxsxs - y{H^XDXs - XsXd^H) + h.c., (14) 


where the contractions of Lorentz and SU{2)s indices are implicit (for Lorentz contractions we use the conven¬ 
tion of Ref. [IS]). Here and henceforth we assnme that xd and Xd couple with equal but opposite strength to 
the Higgs boson,^ and choose y to be real. This prevents contributions to the electromagnetic dipole moment 
of the electron |20|. Furthermore, the interactions in Eq. (141 feature a global SU{2)]i symmetry (broken by 


the hypercharge), which protects the electroweak T parameter from large corrections. After EWSB, the mass 
term for the neutral states is given by Cm = —^{■^o)ijXiXj + h-C-- In the basis Xi = {xs,XdjXz>}> mass 
matrix reads 


(15) 



which is diagonalized by the following transformation 

cos (^4 ®in O 4 

I 

W cos 04 



1 

Am 

= 

Ui) 

\ 


^ Sin 04 

i 

-W cos 04 



(16) 


^ During the final stages of our work, a phenomenological study of the Majorana singlet-doublet model with an analysis of collider 
constraints became available in Ref. |24| . 

^ One could write the Yukawa interaction more generally as Pl^^XdXS ~ VrXsXd^^■ choosing the phase of xs ici Eq. in. 

ms can be made real and positive. Furthermore, by adjusting the phases of xd ^-nd we can make mu and one of the Yukawa 
couplings, pL or i//?, positive and real. Here we assume that p^ = Pr = y, a choice motivated by reducing the corrections to the 
T parameter. Notice that changing the sign of m^) is equivalent to changing the sign of xd (or X%)-> which implies changing 
the sign of pL (or y/?). This results in replacing —m^) with m£) and replacing pvlx/2 with —yu/V^ (or replacing —pvjy/^ with 
pvixpi) in Eq. 1^. 
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The mixing angle 04 is given by Eq. (101 with a 


Xn 


are 


4. The masses of the three neutral eigenstates, ; and 


m 


h,l 


= + TOS ± ATO 4 ), m^ = mD, 


(17) 


respectively, while the mass of the charged state xt) is = rriD as in the Dirac xs case. The couplings of 
the neutral fermions to Z and h are given by 


C D (sin 04 xT “ cos 04 + h.c. 


2c\y 

y 

2 L 


(18) 


sin(204)(x^X° - - 2cos(204)x^X? 


h + h.c.. 


Let us consider the parameter limits in this model. For mu « ms (i.e. \mD — rnsl 'C 2yv, but mo ^ 2j/u), 
the mixing angle is 04 ~ 7r/4. We have three Majorana states split by \yv\, that is ~ mo ± \yv\ and 
= mz). The lightest state, x?) is the DM candidate. All Z couplings to the DM field are thus off-diagonal. 
The DM direct detection signal due to tree-level Z exchange is therefore kinematically forbidden, as long as 
the mass splitting Am 4 is larger than several hundred keV. This requirement is fulfilled for all the benchmarks 
that we consider. At the same time, in the above limit of 04 ~ 7r/4 the coupling of the Higgs boson to the DM 
field X? is maximal and dominates spin-independent DM interactions with nuclei. 

In the limit |mD — ms| ^ yv with mo 3> ms, we have 04 ~ 7r/2. The states x° and Xm are quasi-degenerate 
with masses m° = mo + {yvY/mo and m^ = mo, respectively, and form a pseudo-Dirac state. Its couplings 
to the Z boson are that of the neutral component in a Dirac fermion electroweak doublet, if the mass splitting 
can be ignored. The lightest state x? is mostly a singlet with mass m° = ms — {yv)'^/mo- In the decoupling 
scenario mo ^ ms,yv, both the Z and Higgs exchange in direct DM searches are thus absent at the tree 
level. This feature protects the decoupling scenario from being excluded by direct DM searches. If instead the 
singlet decouples, i.e. if ms 3> mo, then the mixing angle is 04 ~ 0. Now x? and Xm are quasi-degenerate 
states with mass m° « m^ = mo- They form a pseudo-Dirac fermion with unsuppressed couplings to the Z 
boson. This limit is therefore excluded by direct DM searches. 

For mjj < 2yv, the lightest neutral state can be a pure doublet, Xm- In this scenario, DM interactions 
with nuclei are absent at tree level. Up to radiative corrections t he D M state Xm is mass-degenerate with the 
charged state Xo, which leads to strong co-annihilation (see Sec. HI). 


C. The Doublet-Triplet Model 


In this model we assume that the dark sector consists of two fermions, an electroweak doublet, xd, and a 
triplet, XT, that transform under SU{2)z x U{1)y as 


XD ~ (2,1/2-tr), XT ^ {3,r). 


(19) 


We consider the cases r = 0, — 1, in which both xd and xt have neutral components. For r = 0, xt can be 
either a Dirac fermion or a Majorana fermion, while xd is always a Dirac fermion. We thus consider three 
cases: r = — 1 with Dirac triplet, r = 0 with Dirac triplet, and r = 0 with Majorana triplet. 

a. Dirac triplet fermion, r = —1. In this case the electroweak triplet is composed of a neutral state, Xt) 
a singly charged state, Xt? ^cid a doubly charged state, Xt~’ while the electroweak doublet contains a neutral 
state, Xz), and a charged state, Xd^ 


XD = 



( Xt/x^^ Xt ^ 
y Xt ~XtI'^) 


The mass terms and Yukawa interactions are 


( 20 ) 


A -mnXDXD - mr Tr(xTXT) - {vXdXtH -k h.c.). 


( 21 ) 
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After EWSB, the corresponding mass matrices in the bases {Xt^Xd} {Xt’Xd} 


Mo = 





( 22 ) 


respectively, while the mass of Xt ^ 7 ’- EWSB thus introduces mixing between both neutral and charged 
states, so that the mass eigenstates are given by 

x)! = cos6»2Xt+ sin6»2Xi3, X(!) = cos 6»i xj + sin 6»i x^ > ( 23 ) 

X° = -sin6»2XT + cos6»2Xi), xt = -sinflixj +cos6»iXc, 

with the corresponding mass eigenvalues 

'm1i = ^(mD+mT±Am2), = |(mD + my ± Ami). (24) 

The mixing angles 6 * 2,1 and the mass splittings Am 2 ,i are defined in (ITM, with ms —t mx- Since the splitting 
in the neutral sector is larger than in the charged sector, Am 2 > Ami, the lightest neutral state x? is a 
potential DM candidate if mjj and m^ have the same sign. As in Sec. |IIB[ we use the freedom in the phase 
of xd to make y real and positive. In the basis of mass eigenstates, the interactions of the neutral states with 
the Z and Higgs bosons are given by 




1 - 


sin^ 02 \_ 


Cw L 
_ Ji 

V 2 I 


tnl^xl + (1 - 


COS 


-) X?7'^X? 


^(xM + X?7V.) 


(25) 


sin( 202 ) (X^X° - X?X?) + cos( 202 ) {xlx° + X?X°) 


h. 


We consider two parameter limits, i) almost degenerate doublet and triplet, and ii) the decoupling limit. The 
degenerate case occurs for mu ~ mr, so that \mD — mrl ^ | 2 /n|. In this case the mixing is maximal in both 
the neutral and charged sectors, 02,1 — The mass eigenvalues are m° ; ~ mx ± \yv\/^/2 for the neutral 
states, and m((; ~ mx ± \yv\/2 for the charged states. 

In the decoupling limit mo, m^ ^ yv, DM is either predominantly the neutral component of the doublet (for 
mu < mx) or of the triplet (for mx < mu), with a mixing angle of 0[y^v^ j{mx — mu)^)- The masses for the 
neutral states in the two cases are m° ~ "m-xi^x) + ?/^'c^/( 2 |mu — mx\) and m° ~ 'rnx(D) ~ J/^'c^/( 2 |mu — mxl), 
while for the charged states the mass deviation from mD,x is half as large. Since in this model the coupling of 
DM to the Z boson is not suppressed, the direct DM detection searches exclude the entire region of parameter 
space that is consistent with the observed DM relic density. 


b. Dirac triplet fermion, r = 0. The electroweak doublet is composed of a neutral state, Xuj ^ charged 
state, Xuj while the electroweak triplet is composed of one neutral state, Xtj i*wo charged states, Xt^ x'x j 


Xd = 



( XVV 2 


xj 


Xt I /_ 0 / /o 

\ Xt -XT/vSy 


(26) 


The mass terms and Yukawa couplings are given in Eq. (21). After EWSB, the mass matrices for neutral 
and charged states are similar to Eq. ( 22 ), but with Mo ^ M+ and M- —t Mo, in the bases {Xt>Xu} 


{XtiXu}) respectively, while Xx ^ mass of mx- The heavy and light states in the neutral and charged 
sectors have masses 


= 5 (wu + mx ± Ami), = 4 (mu + mx ± Am 2 ). 


(27) 


The lightest neutral state can thus be lighter than the lightest charged state, if mu and mx have opposite signs. 
However, as for the case r = — 1, this model is excluded by constraints from the DM relic density and direct 
detection searches. In particular, obtaining the correct thermal relic density requires sizeable doublet-triplet 
mixing (see Sec. HI), which leads to large direct detection rates from coupling of the doublet component to 
the Z boson. 
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c. Majora na triplet fermion, r = 0. Taking instead the electroweak triplet to be a Majorana fermion, i.e. 
= x!t (26), while the doublet is still a Dirac fermion, the mass terms and Yukawa interactions in the 


two-component notation are (see also Eq. (HI) 


Cm 3 niDXD^XD - ■^mTTT{xTXT) - y{H^XTXD - x'd <^XtH) -f h.c.. 


(28) 


where, as in the Majorana singlet-doublet model, we have introduced two left-handed Weyl fermions trans¬ 
forming under SU{2)l x U{1)y as xd ~ (2,1/2) and Xd (2,—1/2)- As for the Majorana singlet-doublet 

model, the form of the Yukawa coupling terms is restricted by a global S'C/(2)/j symmetry to protect the 
electroweak T parameter. We use the freedom of field redehnitions to make mo and y real and positive.^ 
The mass rriT is in general a free complex parameter, but restricted to real values in our analysis. After 
EWSB, the Yukawa interactions mix the triplet and the doublet components. The three neutral Weyl fermions 
Xi = {Xt> Xc j X%} have a Majorana mass term Cm C> —l(7Wo)ijX/Xj +h.c.. The charged fermions have Dirac 
masses Cm C> —Xmi^ + )mnXt + h.c., where the negatively charged Weyl fermions are Xm = {Xt^x'd}’ 
the positively charged ones, xt = {Xt^Xz)}- The two mass matrices are 


WT 2/2 ~y2 

Mo = I 2 /f 0 -rriD 

-yl -mo 0 


M+ = 


rriT y 




\yT 2 


(29) 


The mass matrix for neutral states is diagonalized by, cf. Eq. (16), 


cos 02 sin 02 sin O 2 


(A 


Xb 

= 

\xl) 

\ 


vA 


vA 



sin 02 cos 02 ^ cos 02 
whereas the charged mass eigenstates are 

xt = cos02xj -bsin02x/), xt = - sin 02xJ + cos 02x/) • 

The mass spectrum is given by 

”T-a.c = '^t,c = 5 (niD + m-T ± ATO2), ml = rriD- 


(30) 


(31) 


(32) 


In the basis of mass eigenstates, the couplings of neutral fermions to the Z and Higgs bosons are given by 


CD -i- 
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2 cw 


2^2 L 


(sin02 Xa* - COS02 Xc*)o-'"Xb^/x + h.C. 
sin(202) (x°Xa - XcXc) - 2 cos(202) x°Xc 


(33) 


h + h.c.. 


This scenario corresponds to the wino-higgsino system (with decoupled bino) in the MSSM for tan /3 = 1 and 
y = g. The DM phenomenology of the doublet-triplet Majorana model for the case where 0 < mu < 200 GeV 
and 2/^1 has been studied in Ref. |26j.'^ In this parameter region, one typically has the mass ordering 
mb < \ma,c\, so that the DM candidate x? = Xb is a pure doublet fermion. In this case x? has no diagonal 
tree-level couplings to the Z and Higgs bosons. The direct detection cross-section in this scenario is loop- 
induced and thus suppressed (see Sec. IV). 


® The sign of y can be flipped by adjusting the phases of xo and without affecting other terms in the Lagrangian. Physical 
observables are thus insensitive to y -H- —y. 

^ Our notation corresponds to the one used in Ref. [26] for Xi = ~'^Xb^ X 2 ~ Xcj X 3 = Xa? Xi” = xi■> X 2 ~ ~Xa , and 
Md = -mo- 











Xi 


Xi 




V\/x/\/' w 


Xi 


Xi 




Xmji 


A/\/\/\/' Z 


%/ 




Xm,h 



FIG. 1: Feynman diagrams for the dominant annihilation channels in the Majorana DM models with fermion mediators. 


On the other hand, if mu, ttit 3> 200 GeV, the mass ordering is rric < rrif, < nia and the lightest neutral and 
charged states are mass-degenerate at tree-level. This degeneracy is lifted by one-loop corrections involving 
gauge bosons, leading to [27U^ 


M+ = 


rriT 


6 m^ 


'V2 


yT2 

ttid + Sm'jj 


(34) 


In the limit mT,mD, \mT — mz the radiative splittings are given by 


Smt = —{mw - c^mz), 

OTT 


OTT 


(35) 


The corrections to the off-diagonal elements in the mass matrix have been neglected above, which is justified 
for yv <C \mT — mul- The one-loop corrections can also be neglected for the calculation of the mixing angles. 
Since the corrections Sm^ and Srn^ are positive, the lightest state of the spectrum is the neutral state Xc> ^ 


DM candidate. As is apparent in Eq. (33), the Z boson coupling to a DM pair jqxZc is absent and the Higgs 


coupling is proportional to sin( 202 )- Direct detection signals are thus suppressed for a small mixing angle 02 - 


III. THERMAL RELIC DENSITY 


As described above, we assume that the fermionic Higgs portal is responsible for explaining the entire dark 
matter density through thermal freeze-out in the early universe. 

In the singlet-singlet model, DM preferentially annihilates into ^ 2/12 and h 2 hi final states, if kinematically 
allowed. The amplitude for XX ^ 2^2 is proportional to cos^ a, where a is the h-S mixing angle, so that 
annihilation can be efficient even for very small values of a. For < (m/^ -|- 771 ^, 2 ) 72 , the main annihilation 
channels are into hihi and W~^W~, where the latter proceeds via off-shell hi 2 exchange in the s-channel. For 
small values of the bb final state can also become relevant. The rates for these processes grow with sin a. 
The requirement of a sufficiently large annihilation cross-section then imposes a lower bound on sin a. 


For the Majorana DM models with fermi on me diators (the Majorana singlet-doublet model in Sec. HB and 
the Majorana doublet-triplet model in Sec. HC), the main channels for pair annihilation of the neutral DM 


candidates involve WW and ZZ final states. These processes are mediated by one of the mediator fermion 
states in the t-channel, see Fig.fH since the Zx^X^ coupling vanishes exactly. For the Majorana singlet-doublet 
model, annihilation via the s-channel Higgs-boson resonance is also a viable option for « rnhj'l. In this 
case, resonant enhancement from the Higgs-boson propagator leads to a sufficiently large annihilation cross- 
section to produce the correct relic density. The dominant annihilation final states are then given by the 
leading Higgs decay modes, i.e. bb, WW*, gg and r+T“. 

In contrast, in the Dirac singlet-doublet model, the DM annihilation mainly proceeds through s-channel 
Z-boson exchange. Only for very large DM masses, ^ 0(1 TeV), annihilation into WW and ZZ final 
states through t-channel fermion exchange becomes important. 

In the singlet-doublet models (both for the Majorana and Dirac cases), the lightest neutral fermion is 
constrained to be mostly singlet, to avoid the strong direct detection bounds for doublet dark matter (see 
Sec. IV). However, the singlet nature of DM in these models also suppresses the annihilation cross-section. 


thus typically yielding too large of a relic density. Nevertheless, the correct DM density could still be obtained 
if and X^Xm CO- annihilation processes contribute at a sizeable level. As a result, the allowed parameter 

space is limited to relatively small values for the mass splitting mu — ms- One exception is the Higgs resonance 
region for the Majorana singlet-doublet model, where the correct value for the annihilation cross-section can 
be obtained without co-annihilation. 
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FIG. 2: Range of Lagrangian parameters mo and my compatible with the thermal relic density constraints from 
Eq. (361 in the Majorana doublet-triplet model. The curves correspond to different values of the Yukawa coupling y. 


In the scenario with pure doublet Majorana DM, Xmi the degeneracy of the states Xm and leads to strong 
co-annihilation. The correct thermal relic density is obtained only if all dark particles lie above the TeV scale. 
Since such a scenario is not accessible at the LHC or planned future colliders, e.g., a proton-proton collider 
with y/s < 100 TeV or an e+e” collider with ^/s < 1 TeV, we do not investigate this possibility further. 

We have computed the relic density using MicrOMEGAs 3.6.9.2 [31], which automatically incorporates co¬ 
annihilation processes and three-body final states through off-shell W/Z production. The models described in 
the previous section have been implemented into model hies for CalcHEP [32j . which is used for the matrix- 
element generation within MicrOMEGAs. The relic density is then required to match the value obtained by the 
Planck collaboration |33|, 


= 0.1199 ± 0.0022. 


(36) 


For the models with fermion mediators, this reduces the three-dimensional parameter space, e.g. mu, ws/Ti 
for each model down to two independent parameters. This is illustrated by the solid colored lines in Figs. |5}|^ 
It is worth pointing out that the Majorana doublet-triplet model has two distinct regions of parameter 
space that are compatible with the relic density constraint from Eq. (36). The hrst region is realized for 
rriD < 200 GeV and y > 1, leading to a DM candidate = Xb that is an unmixed doublet state with mass of 
a few 100 GeV [2H] . 

The second possibility is obtained for tod, mr ^ 1 TeV. In this case, the lightest ne utral and charged states. 


Xi = Xc Xi = Xc j spht only by small radiative corrections (see the end of Sec. IIG). Asa result, there 


is strong co-annihilation and the relic density comes out too small if the DM mass is below about 1 TeV. 

Fig. 15 shows the parameter space that is compatible with the relic density constraint in both cases, as a 
function of the mass parameters mD,T and for several sample values of y. To ensure perturbativity, an upper 
limit ?/ < 3 has been imposed. Unfortunately, the second scenario, where all the new fermion states are beyond 
1 TeV, is inaccessible at the LHC and planned future accelerators. Therefore we will not study it any further 
in Sec. [VD 


IV. DIRECT DETECTION 

The non-observation of DM scattering off atomic nuclei, referred to as direct detection, leads to strong 
bounds on our models of fermion DM. At the tree level, spin-independent interactions of the lightest neutral 
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Model 

DM state 

G| 

Gl 

fn/fp ~ 

Majorana xs + S 

XS 

- 


1 

Dirac xs + Xd 

XS 

9AT--2si^Qq) 2 n 


1 

cos 

1 

1 

Majorana xs + Xd 

XS 

- 

4'm^ ^ V 

1 

Majorana Xt + Xd 

Xd 

- 

Q + rk^—Sy 

4m“ my/ ^ 

1 

Majorana Xt + Xd 

Xt 

- 

4m^ myy y/2 '' ' 

1 


TABLE I: Effective DM-quark interactions from Z exchange (G%') and Higgs exchange (G^). The column “DM state” 
indicates the dominant component of the lightest neutral state in each model. The ratio fn/fp gives the amount of 
isospin violation, assuming that Z exchange dominates. Furthermore, rrih is the Higgs mass, g is the weak coupling 
constant, and and Qq are the weak isospin three-component and the electric charge of the quark q. The one-loop 
correction 5y is given in Appendix [X] 


state of the spectrum, with quarks q inside the nucleus are generally mediated by the Higgs boson and, 
in case Xi is a Dirac fermion, also by the Z boson. Since the momentum transfer in DM-nucleus scattering 
is much smaller than the Higgs or Z mass, spin-independent DM-quark interactions can be described by an 
effective Lagrangian, 


^EFT D -G| (x? 7 mX°) (97'"'?) - Gl (X°X°) (ot) ■ 


(37) 


The Z- and Higgs-mediated effective couplings and G^ are listed for our models in Tab. We show 
only models with a suppressed Z coupling, which can pass the limits from direct detection searches and 
simultaneously provide the correct DM relic density. By comparing the dependence on the fermion mixing 
angle 9a with the definition in (101, it is apparent that all models favor a small Yukawa coupling y to evade 


bounds from direct detection. At zero momentum transfer, the cross-section for spin-independent scattering 
of a DM particle off a nucleus with mass number A and proton number Z is given by 


a A = 




Zfp + iA- Z)fp 


12 fp=U 


a. 


(4 ,2 


(38) 


where fp and /„ are the DM couplings to protons and neutrons, Op = ky,pfp/'!r is the DM-proton cross-section, 
Mi = rnfrnf /(rn^ + rrii)^ is the reduced mass, and k = 1(4) if DM is a Dirac (Majorana) fermion. In terms of 


the effective interactions in Eq. (37), the DM-proton coupling is given by 


fp = S 


<? 


DZ G% 


E/: 


(p)^g 'ffp 




nip 


mg 


(39) 


and analogously for /„ with p ^ n. Here Sd = 1 (0) for a Dirac (Majorana) DM fermion, q = u,d,s and 
Q = c,b,t denote light and heavy quarks, and = 1 — d s frq^ where nipf^^ = {p\mq qq\p) describes 
the matrix element of quarks inside the proton. 

The results of direct detection searches are usually expressed in terms of the DM-nucleon cross-section 
apf = kp’jqffj/'K^ taking the isotope abundances in the detector material into account and assuming isospin- 
conserving couplings fN = fp = fn and = /i^ « However, if DM is a Dirac fermion, the interaction 
through Z exchange can induce significant isospin violation, resulting in fp ^ /„. When deriving bounds on 
the DM-proton cross-section ap in our Dirac DM models, we therefore rescale the experimental bounds on cr at 
from Ref. [34j 


J2^'n^9A^Z + {Ai- Z)fn/fp\‘^' 


where pi denotes the abundance of isotope Ai. The ratio fn/fp for our models is listed in the last column 
on Tab. |l] The Z-exchange contribution in the Dirac singlet-doublet model is so large that it dominates the 
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DM-quark interaction. The resulting isospin breaking leads to a cancellation in the DM-nucleus scattering 
amplitude. The bounds from direct detection on this model are therefore much weaker than they would be 
when neglecting isospin effects. Since isospin violation in Higgs exchange is very small, the cross-sections in 
all other models are not subject to isospin breaking effects. 

In our numerical analysis, we compare the DM-nucleon cross-sections computed with MicrOMEGAs 
3.6 .9.2 m with the latest bounds from the LUX experiment |3S] and projected bounds for the XENONIT 
experiment [Ml- We use the values for the scalar nucleon quark form factors given in Tab. 3 of Ref. [5T|^. 

For the Dirac singlet-doublet model we take into account isospin violation using (40), with the isotope abun¬ 
dances for Xenon listed in Tab. II of Ref. [51]. The resulting bounds on the parameter space in each model 


are displayed in the figures in Sec. VI 


In models with Majorana DM, Z-mediated spin-independent scattering off nuclei is generally absent at 
tree level, since the Majorana fermion is its own anti-particle, while the vector current is odd under charge 
conjugation. In the Majorana doublet-triplet model with pure doublet DM, also the Higgs-mediated 

scattering vanishes at tree level. The reason is that yJJ does not obtain its mass from the Higgs mechanism and 
does not mix with the triplet through Yukawa interactions. A DM-Higgs interaction is induced at one-loop 
through the exchange of dark fermions and electroweak gauge bosons. 


-y XbX°/i + h.c.. 


(41) 


The Yukawa correction Sy is a function of the parameters y, , and . An analytic expression 

is given in Appendix [A] This vertex generates an effective DM-nucleon scalar interaction through Higgs 
exchange, as given in Tab. |T] Scalar interactions are also induced by box diagrams with W and Z bosons. 
These contributions are numerically sub-leading if the Yukawa coupling y is large, which is the case if DM is 
a doublet. We therefore neglect them in our analysis. 


V. HIGGS DECAYS 

In the singlet-singlet model, mixing between the doublet Higgs h and the singlet scalar S leads to a universal 
suppression of all couplings of the SM-like scalar hi, 

Kt = Kb = Kr = Kw = Kz = COStt, (42) 

where Ki = ghiu/gf^- Thus the value of the mixing angle a can be constrained from measurements of the 
Higgs production rates at the LHC and future colliders iiniiii]. 

In the doublet-triplet model, the decay rate of the Higgs boson into two photons, r(/i — >• 77), is changed at 
the one-loop level due to the presence of virtual charged fermions. The decay ratio with respect to the SM 
rate, Tsuih^ 77 ), is given by 


T{h 77 ) 


rsM(/i^77) 


1 A ^ 
1 -f ^ 


A: 


SM 


Here the one-loop amplitudes for the SM, Asm, and the charged fermions, A^, are defined as 


Asm = ^ NcQfAp(Tf) + Ab(tw) and A^ = '^ Q^y^^Ap(T^) 


(43) 


(44) 


5 


It should be kept in mind, however, that the values for are subject to large theoretical uncertainties |371I38| . Hadronic 

uncertainties on and can be reduced by using the method described in Ref. m- 
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where Ti = rn\/(Aml) and the loop functions for r < 1 are, 

2 
T 


— i arcsin 


^b{t) = —^| 2 t^ + 3 t + 3(2r — 1) arcsin^ ^/t 


(45) 


In terms of the mass eigenstates from Eq. (32), the amplitude is given by 


— 


j/sin(26»2) r V 


V2 


{: 


1T 


^p’(t^+)}, sin(26i2) = 


V2yt 


mf - mt 


(46) 


Notice that is proportional to and thus independent of the sign of the Yukawa coupling. Furthermore, 
the contribution vanishes in the limit of no fermion mixing, so that at least two charged fermions are required 
to induce the Higgs decay. Since Ap(t) is monotonously increasing with r and since per definition, 

is always negative for positive The contribution of charged fermions thus depletes the —)■ 77 

the measurements of Rj by the LHC 


decay rate below the SM expectation. As we will discuss in Sec. VI 


collaborations set strong constraints on the Majorana triplet-doublet model in the parameter space where the 
Yukawa coupling is sizeable. 


VI. RESULTS AND COLLIDER BOUNDS 


In what follows, numerical results for the allowed parameters for each of the Higgs portal models introduced 
in Sec. [H] will be shown. Besides constraints from the relic density and direct DM detection, we also evaluate 
bounds from current LHC data at ^/s = 8 TeV (LHC 8 ), as well as projections for upcoming LHC runs 
at v^s = 14 TeV (LHC14) and future planned colliders. In particular we consider the reach of a proton- 
proton collider with y/s = 100 TeV (called FCC-hh in the following) and of a high-energy e+e” collider with 
^/s = 1 TeV, such as the planned International Linear Collider (ILC). 

In the singlet-singlet model, the experimental sensitivity to the new states is due to the Higgs-singlet mixing. 
This modifies the couplings of the SM-like Higgs, hi, to the SM fermions and gauge bosons, which can be 
determined by measurements of the Higgs production rate in different decay channels. One can also search 
directly for the second scalar state, / 12 , in the same channels as for the SM Higgs, but with a different mass 
and different production cross-section. The most sensitive channels are /12 WW and /i 2 ZZ. 

In the singlet-doublet and doublet-triplet models, the particle spectrum involving several neutral and charged 
fermions can lead to LHC signatures with leptons and missing energy. The most important channel is given 

by 


qq' —^ W* —>■ X Xm.ft) 

X X?^ Xm,h ^ Xi Z* ^ (^ = e,/r) 


(47) 


and its charge-conjugated versions. The production process, mediated by an off-shell W boson, is present both 
in the Majorana and Dirac DM models, and its rate is governed by the SU(2) gauge quantum numbers of the 
fermion multiplets. The only decay channel for the lightest charged fermion is also through an off-shell W 
boson. For the next-to-lightest neutral fermion, tree-level decays can be mediated by an off-shell Z or Higgs 
boson, but the latter is suppressed to a negligible level by the small Higgs-boson width. 

The process (471 leads to a final state of three leptons (3^) and missing energy, with a phenomenology similar 
to the neutralino-chargino sector of the MSSM with heavy scalars. This signature has been searched for by 
the ATLAS and CMS experiments [HI |33] . If the mass difference between the DM fermion and the heavier 
fermions is small, which is the case in the regions of parameter space where co-annihilation is important, the 
leptons can become too soft to pass the usual multi-lepton selection cuts. In this case, mono-jet searches 
provide the strongest limits from the published 8 -TeV LHC data |44l [45] . 


Notice that the heavier neutral fermions Xm h general also decay through a loop-induced radiative 


decay, x° 


X° 


7 . However, the branching fraction is suppressed to the percent level if either the heavy 


state, Xm hy ^be light state, X/j has a dominant doublet component 
considered here. 


as is the case for the models 
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A. Singlet-Singlet Model 


At colliders, this model can be most effectively tested by searching for the effects of the singlet scalar 
mediator, S. It mixes with the Higgs, h, so that the mass eigenstates hi and /i 2 are admixtures of S and h, 
see Eq. ([^. As a result, the production and decays of the SM-like scalar hi to SM fermions and gauge bosons 
are suppressed by a common factor cos^ a, as discussed in Sec. 

When interpreting the CMS and ATLAS Higgs measurements, we need to distinguish two cases. If < 
mhj^l2, then hi XX is allowed, and thus the invisible branching ratio of the Higgs, Br(h —inv), could be 
sizeable. Allowing for an invisible decay of the Higgs and demanding Ky < I, ATLAS obtains in a global 
analysis Ky > 0.93 at 95% C.L., Kf = 1.05 ± 0.16 and Br(/ii —inv) < 0.13 |47j. (CMS obtains a much looser 
bound Br(/ii —>■ inv) < 0.49 for Ky < 1 [IS].) The invisible decay width is given by 


r(hi ^ xx) 


"m-hi 2 







4to^\ 1/2 


(48) 


and thus depends on Sa, as well as on y and j/ 5 . The correct relic abundance can be obtained for y as small as 
0(10“^) [TS]. This means that, allowing for all values of y, j /5 that satisfy the Planck measurement (36), the 
strongest constraint on Cq, is not the bound on Br(/ii —>■ inv), but rather the bound on ny- We therefore have 
at 95% C.L. 


Ca > 0.93 for < mhj^j2, 


(49) 


or Sa < 0.37. 

If > mhil2, the invisible decay width of the Higgs is zero. We can therefore use the results of global 
fits of the Higgs data, where only Kf = Kb = Kr = nt and Ky = kw = kz are varied. Combining the CMS 
results, Ky = l.OI ± 0.07 and k/ = 0.87lo;^3 [15], and ATLAS results, Ky = 1.09 ± 0.07, Kf = l.ll/to^is HZ], 
gives Cq = 1.035 ± 0.045 or a 95% C.L. lower bound 


Cq > 0.948 for 


(50) 


corresponding to Sq < 0.318. The bound (50) also applies for just slightly below m,h^l2, where the invisible 
decay width of the Higgs is negligible because of phase space suppression. The above results assume that /12 
is not (approximately) mass degenerate with hi. In the remainder of this subsection we will assume that this 
is true and that, furthermore, ^2 is appreciably heavier than hi. For bounds away from this limit in a subset 
of our parameter space, see for instance Ref. |49) . 

The bounds on a can be improved by Higgs coupling measurements at future colliders |41j . The LHC with 
i/s = 14 TeV and 3000 fb“^ luminosity (HL-LHC) can establish a limit of Cq > 0.9 at 95% C.L., assuming 
that the central values of all observables agree with the SM prediction. The projected reach of ILC with 
ys < 500 GeV and 500 fb"! is eg > 0.96. 

In Fig. [^ these bounds are compared to the region of parameter space of the singlet-singlet model that 
is consistent with the relic density constraint in Eq. (36). The points in the figure are obtained through a 


parameter scan, where the model parameters are varied randomly within the following ranges: 


100 GeV <my< 1000 GeV, 


- 10 < A' < 10, 


- 1000 GeV <ns < 1000 GeV. (51) 


The values of y' and ms are determined from the mixing angle Sq and the masses mh^ see Eq. ([^, while 
A 5 is irrelevant for the relic density calculation. For simplicity and to avoid constraints from GP violation, 
we take 2/5 = 0. For each random point, the value of y is fixed by the relic density constraint, Eq. (36), and 


the condition j/ < 3 is imposed to ensure perturbativity. Points shown as yellow circles are excluded by limits 
on the spin-independent direct detection cross-section from LUX [35], while the red crosses and blue dots are 
within and beyond the projected reach of XENONIT |5B], respectively. 

One can also search directly for the heavy Higgs, / 12 . The decay widths to SM fermions and gauge bosons 
are given by 


r(h2 ^ //) = slr{H ^ //) 


SMj 


r(/i2 ^ vv) = slr{H ^ vv) 


SM, 


(52) 


where V = W,Z, 7 , g, and r(i7 —)■ ArA')sM is the partial width of the would-be SM Higgs if it had a mass mh 2 ■ 
If only the SM decay channels are open, the branching ratios of /12 are thus given by the branching ratios of 
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FIG. 3: Allowed parameter space for the singlet-singlet model consistent with the requirement that the thermal relic 
density of x accounts for all dark matter in the universe, — Udm- The different colors (shapes) of the points indicate 
current and future 90% C.L. exclusions from direct detection experiments. Also shown are current 95% C.L. limits from 
Higgs coupling measurements at the LHC (dashed line), and future projections for LHC 14 with 3000 fb“^ (dotted) and 
ILC with ^/s < 500 GeV (long dashed). The left panel corresponds to < mh.il'i-i which forbids the annihilation 
channels XX ^ 2 / 11 , 2 , while in the right panel > mh^l2. 
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the SM Higgs with mass rrih^. Similarly, the production cross-section is < 7 ( 112 ) = s'^a{H)sM, and is given by 
the would-be SM Higgs production with mass . 

For mh 2 > decay /12 —>■ hihi is kinematically allowed. It proceeds through the interactions 


-Sh^ - X'vS^h - fisS^ 

V 2 


= h2hj 


3 — Sacl - ^{cl- 2sl)ca - X'v{sl - 2c^)sq - 3fJ.SCasl 

V 2 




(53) 


If the mixing angle is small, Sa 1, then the / 12 / 11/11 coupling is equal to n'/2. This means that the /12 hihi 
branching ratio can be large and will dominate over /i 2 —f H^H^, ZZ, ti for fj! ^ v. On the other hand, if /r' ^ v, 
the decay /12 —>■ hihi will be subleading. It is also possible to make the /12 —>■ hihi branching ratio small while 
keeping Sa large, by canceling different contributions in Eq. (53). For rrih^ > 2.m^, the decay /12 —t inv. is open 
and can become dominant for sizeable values of y oi y^. The actual constraint from /12 decays thus strongly 
depends on the specific realization of the model. 

In Fig. 1^ (left), we show the bounds from searches for direct decays of /12 into SM final states. They 
can be expressed as a constraint on x Br(h 2 —>■ VV,ff) as a function of where V = W,Z,^,g. 
The most constraining channels are ^2 —f WW, ZZ, with the resulting bound from CMS shown as a dashed 
black curve [52]. The searches in the di-photon channel, /12 —t 77, are effective at lower /12 masses. The 
resulting bounds from ATLAS [50| and CMS |5T| are shown as solid curves in the upper left corner. In 
Fig. (right), we also show bounds on x Br(/i 2 —>■ hi hi) from di-Higgs searches. The CMS [S3] and 
ATLAS |S3] bounds from h 2 —>■ hihi —)■ 2^2b are shown as solid black and blue lines, respectively. The CMS 
bound from h 2 —t hi hi —)■ ib |55| is shown as a dashed black line. There are also CMS di-Higgs searches 
with di-photon and leptonic final states m and h 2 —>■ tt |58j . However, with their current precision these 
results do not constrain Sa- We see that relatively large mixing angles Sa ~ 0.3 are allowed for all h 2 masses, 
comparable to the constraints from global Higgs coupling measurements in Eq. (49). In Fig. 13 we also show 
as dashed (solid) red lines the projected exclusions at 14 TeV LHC with 300 flv^ (3000 fb'^) obtained in 
Ref. [S5]. At the end of the high luminosity LHC run, mixing angles as small as ~ 0.05 can be probed for 
mh 2 ^ 400 GeV. Notice that for large h 2 masses, large mixings are additionally constrained by electroweak 
precision tests HH EH] . 

At a 100 TeV collider, one can also search directly for DM production through an off-shell singlet mediator 
in the monojet signal, pp h*(—>• xx)j, even if hi does not decay to DM [10]. This signature is quite 

challenging due to its very small cross-section, and can thus be observed only in a small parameter region with 
just above mh^/2. More promising is the signal with two jets and missing energy, where the presence of 
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FIG. 4: Left: Existing bounds on s% in the singlet-singlet model from h 2 —>■ 77 (solid lines from ATLAS |50| and 
CMS [^) and from h 2 —> WW, ZZ (dashed black line, CMS [^) as a function of the heavy scalar mass, Right: 
Bounds on x Br(/i 2 —>■ hihi) from /12 —> hihi —> 2^2b (solid black from CMS |53| . solid blue from ATLAS [54]) 
and from /12 hi hi —>■ 4b (dashed black from CMS [55]). The projected exclusions at the 14-TeV LHC with 300 fb“^ 
(3000 fb“^) [^ are shown as dashed (solid) red lines. The colored points indicate the parameter region consistent 
with the relic density constraint, = Udm, with the different colors (shapes) denoting current and future 90% C.L. 
exclusions from direct detection experiments. 


an additional jet allows for easier discrimination from the background |60j . 

As can be seen from the figures, current and future direct detection experiments can cover a significantly 
larger portion of the parameter space of the singlet-singlet model than collider experiments. However, there 
are also many parameter points that are not constrained by available data from LUX, but that are excluded by 
LHC results for Higgs couplings and heavy Higgs searches. Therefore, collider and direct detection experiments 
provide complementary information for probing this model. On the other hand, the correct relic density can 
be realized for very small values of the mixing angle, sin a ^ 0.01, which cannot be tested conclusively either 
at colliders or through direct detection. 


B. Majorana Singlet-Doublet Model 


For > 100 GeV, co-annihilation is required in the Majorana singlet-doublet model to obtain the correct 
thermal relic density for weak Yukawa couplings, y < 1. This implies that the mass difference is 

preferred to be a few tens of GeV. For DM masses in the TeV range the model thus becomes significantly 
fine-tuned. We therefore restrict ourselves to the range m° < 1 TeV. 

If the mass difference becomes too small, the relic density constraint cannot be satisfied even for very small 
values of y. This is depicted by the shaded region at the bottom of Fig. On the other hand, large values of 
y are constrained by limits on the spin-independent direct detection cross-section from LUX [3S| (see shaded 
region at the top of the plot). 

In the remaining allowed part of parameter space, the mass difference — rrif is relatively small, so that 
the leptons from the process (47) are soft and fail the selection cuts for the 3£ signature [321 03] . On the 


other hand, the production cross-section for (471 is too small to be constrained by the available mono-jet data 
[44l l45] . As a consequence, no bound on the cosmologically preferred parameter space of the singlet-doublet 
model is obtained from LHG8 data (see also Refs. [3TH55] ). Similarly, the projected reach for the M signal of 
the LHG14 [53105] and FGG-hh [551EZ] does not extend into the white region in Fig. 

Instead, requiring a hard initial-state jet can help to trigger on events with soft leptons and improve the 
signal-to-background ratio for this case. Several authors have analyzed this signature, consisting of at least 
one hard jet, large missing energy, and at least two soft leptons, and found it to be promising for the parameter 
region preferred by co-annihilation [3i[5Il[5BH7n]. We have obtained the estimated 95% G.L. reach of LHG14 
by recasting the analysis of Ref. [61]. Goncretely, Ref. [5T] contains results for a Majorana singlet-triplet 
scenario, with the dominant signal contribution stemming from the process (47). In our case, since we have 
heavy doublet fermions instead of a triplet, the production cross-section is reduced by a factor of two. After 
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FIG. 5: Allowed parameter space for the Majorana singlet-doublet model, under the assumption that Q^o = Udm- 
The solid lines indicate the correct DM density for different values of the Yukawa coupling y. The shaded regions are 
excluded by yielding too small of a relic density (bottom) and by direct detection limits from LUX (top, red dashed). 
Also shown are projected 95% C.L. limits from jet plus soft leptons at LHC14 with 300 fb“^ (short dashed) and 
3000 fb“^ (dotted) and FCC-hh with = 100 TeV and 3000 fb“^ (dot-dashed), as well as from ILC with y/s = 1 TeV 
(long dashed). The red dotted line depicts the expected 90% C.L. direct detection limit from XENONIT. See text for 
details. 


accounting for this change in the production rate, the projected limits shown by the short-dashed (for 300 fb“^) 
and dotted (for 3000 fb“^) lines in Fig. are obtained. As can be seen from the figure, such an analysis will 
only be able to test this model for relatively small DM masses, < 250 GeV, and only in the high-luminosity 
run of the LHC. 

However, a much larger part of the parameter space can be probed with FCC-hh with ^/s = 100 TeV. We 
use the results for higgsino and higgsino-bino scenarios in Ref. |69j . which directly correspond to our Majorana 
singlet-doublet model. The authors of Ref. [69] have simulated the signals and backgrounds for the monojet 
signature and the jet plus soft leptons signature. The projected 95% C.L. limit is shown by the dash-dotted 
line in Fig. We find that the astrophysically favored parameter space for < 900 GeV can be covered. 

Also shown in the figure is the projected reach of the ILC with ^/s = 1 TeV, which has been evaluated 
by extrapolating the results of Ref. m- Due to its clean environment and low backgrounds, the ILC will 
be able to probe this model for masses of the doublet fermions, = m’*' = mjj, up to almost half of the 
center-of-mass energy (see long-dashed line in Fig. [^. 

In addition to constraints from future collider experiments, the allowed parameter space of the Majorana 
singlet-doublet model will also be probed by upcoming direct detection experiments. The dotted red line in 
Fig. El indicates the projected reach of the XENONIT experiment jSSj- As evident from the figure, it will be 
able to cover large parts of the parameter space, except for very small Yukawa couplings, y < 0.02. 


The correct relic density in this model can also be achieved by annihilation of x? pairs through the Higgs 
resonance. In this case, co-annihilation does not play any role, and the other fermion states can be much 
heavier than the DM state, as shown in Fig. In this scenario, both the annihilation and the direct detection 
cross-sections are mediated by Higgs exchange, so that the direct detection bound from LUX [3S] becomes 
independent of the heavy fermion mass in this plot. It excludes part of the parameter space that produces the 
correct relic density. 

The remaining parameter space can be probed in three ways. The M signature from the process (47) can 
be observed if the mass difference — m° is sufficiently large, but the heavy states Xm are not too 

heavy. The first condition ensures that the leptons can pass the trigger and selection requirements, while the 
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FIG. 6: Allowed parameter space for the Majorana singlet-doublet model in the Higgs resonance region, with the 
constraint = Hdm- The solid lines indicate the correct DM density for different values of the Yukawa coupling y. 
The red-shaded region is excluded by direct detection limits from LUX. Also shown are projected 95% C.L. limits from 
31 searches at LHC 14 with 3000 fb“^ (dotted horizontal lines and shading) and FCC-hh with 3000 fb“^ (dot-dashed 
horizontal line and shading). The vertical lines depict the projected limits from h —>■ invisible measurements at LHC 14 
with 300 fb“'^ (short dashed) and 3000 fb“'^ (dotted), and at ILC with ^/s = 250 GeV (long dashed), as well as from 
direct detection searches at XENONIT (red dotted). See text for details. 


second is related to the need for a large enough signal production cross-section. For LHC14 with 3000 fb“^, 
we estimate its 95% C.L. exclusion limit for this signature based on simulation results from the ATLAS 
collaboration [^. To account for the smaller cross-section in the Majorana singlet-doublet model compared 
to the supersymmetric singlet-triplet scenario studied by ATLAS, we conservatively take the 5cr rather than 
the 95% C.L. contour from Fig. 10 in Ref. [65].® For FCC-hh with ^/s = 100 TeV and 3000 fb“^, we adopt the 
results from Sec. III.C in Ref. [66]. As evident from the horizontal bands in Fig. the combination of LHC14 
and FCC-hh results can cover most of the astrophysically allowed parameter space, except for relatively large 
values of the Yukawa coupling y. 

Alternatively, this scenario can be tested by measurements of the invisible Higgs width. For y® < mh/‘2, the 
Higgs boson can decay into pairs of DM particles, h —)■ y®x®, which escape undetected. Future LHC data will 
be able to put strong limits on the invisible branching fraction of 17% with 300 fb“^ and 6% with 3000 fb“^, 
under somewhat favorable assumptions for the systematic errors |64j . A much more precise constraint of 
Br(/i —)■ inv) < 0.3% is expected from ILC data taken at yfs = 250 GeV (see Tab. 2.6 in Ref. [72]). The 
corresponding bounds on the parameter space of the model are depicted by the vertical lines in Fig. [^ These 
bounds do not depend on the heavy fermion mass, since both the thermal relic density and the Higgs 
invisible width are governed by the same hxiX^ coupling. 

Finally, future direct detection searches by XENONIT [36] will probe most of the allowed parameter space, 
as depicted by the red dotted lines in Fig. [^ Combining the projected limits from collider and direct detection 
experiments, only the parameter region where the DM mass is very close to half the Higgs mass, 60.5 GeV < 
TO® < 62 GeV, and the mediator mass is very large, toJ^ > 3 TeV, will remain inaccessible. 


® While this choice will underestimate the reach of LHC14, we note that the covered parameter space does not change much 
between the 95% C.L. and 5cr contours from Fig. 10 in Ref. m- 
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FIG. 7: Allowed parameter space for the Dirac singlet-doublet model, with the constraint = Ddm- The solid 

lines indicate the correct DM density for different values of the relative mass difference 5m = {'mil ~ mil)/ml- The 
red-shaded region is excluded by direct detection limits from LUX. Also shown are projected 95% C.L. limits from 
jet pins soft leptons searches at LHC14 with 300 fb“^ (short dashed) and 3000 fb“^ (dotted), and at FCC-hh with 
3000 fb~^ (dot-dashed). The reach of ILC with -^i = 1 TeV is depicted by the long-dashed line, while the red dotted 
line indicates the expected 90% C.L. limit from XENONIT. The gray region is excluded by the Fermi-LAT indirect 
DM searches. 


C. Dirac Singlet-Doublet Model 

As for the Majorana case, co-annihilation in the Dirac singlet-doublet model is active for most of the 
parameter points that satisfy the relic density constraint. In fact, the predicted value for in this model is 
very sensitive to the mass difference ml — ml. Therefore, instead of displaying our results in the plane of ml 
and ml — ml, we chose ml and the Yukawa coupling y as independent variables in Fig. ItI 

As expected, large regions of the parameter space are excluded by direct detection limits from LUX 
(shaded region in plot), since the doublet component of xl couples to the Z boson and thus leads to sizeable 
DM-nucleon interactions. In the remaining part of the parameter space, the mass difference m° — ml is small, 
which is illustrated by the colored solid curves in Fig. This bound will be improved by future results from 
XENONIT [3S] (red dotted line in plot), but the region with small mass differences and small Yukawa couplings 
will remain difficult to probe by direct detection experiments. The upper left corner of the parameter space 
in Fig. with large y and small ml is also excluded by indirect searches. The 95% C.L. bound from the latest 
analysis of dwarf spheroidal galaxies by the Fermi satellite m is shown as a long-dashed line. The gray region 
to the left of the line is excluded. 

As explained in the previous section, scenarios with such small mass differences can be best searched for by 
using the hard jet plus soft leptons signature at the LHC. The expected 95% C.L. reach of LHC14 with 300 fb“^ 
and 3000 fb“^ from the analysis in Ref. [5T] is shown by the dashed and dotted lines in Fig. respectively. 
These bounds have been obtained by rescaling the results of Ref. m to account for the production cross-section 
of pp —>■ X~Xh^ in Dirac singlet-doublet model. Also shown are the projected reach of FCC-hh for 

the jet plus soft leptons signal from the analysis of Ref. [69] (dash-dotted line), as well as the reach of ILC 
with -y/s = 1 TeV (long-dashed line). 

As is evident from the figure, LHC 14 will be able to probe part of the astrophysically viable parameter 
space, but only for light DM with m° < 250 GeV. In contrast, ILC with ^/s = 1 TeV can extend this reach 
to about ml < 490 GeV, while FCC-hh is expected to cover the entire allowed parameter region. 
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FIG. 8 : Allowed parameter space for the Majorana doublet-triplet model, with the constraint = Udm- Each point 
in the m^-y parameter plane has two solutions for 5m = — m^)/m^, which are shown in the two panels. The 

red-shaded region at the bottom right is excluded by yielding too small of a relic density, while the light shaded regions 
at the top are constrained by data from LUX. Also shown are projected 95% C.L. limits from 3£ searches at LHC14 
with 300 fb“^ (short dashed) and 3000 fb“^ (dotted), and the 90% C.L. reach of XENONIT (red dotted). In the left 
panel, the entire allowed parameter region can be covered by LHC with 3000 fb“'^, as well as by XENONIT. Notice 
that the whole viable parameter space in both panels is already excluded by h —^ 77 data from LHC 8 , assuming that 
there are no new fields coupling to the Higgs boson beyond the doublet and triplet fermions. 


D. Majorana Doublet-Triplet Model 

In this subsection, we explore the parameter region of the Majorana doublet-triplet model for a pure doublet 
DM candidate, corresponding to the left panel of Fig. This scenario is characterized by sizeable mass 
splittings between different dark sector fermion states, so that co-annihilation is not important. 

For this model, it is not possible to display our results in the plane of and — m°, since for some 
values of and there are more than two possible values of y that satisfy the relic density constraint. 
Instead, as in the previous subsection, and the Yukawa coupling y are used as independent variables. Even 
in this case, there is a two-fold ambiguity for the value of tot for each point in the m^-y plane, as can be seen 
in Fig. (left). These two solutions are shown in the two panels in Fig. The solid colored curves indicate 
the correct relic density for different values of the relative mass difference 5m = — m°)/m°. Notice that 

large values of y are constrained by vacuum stability [55]. 

The region with small Yukawa couplings is excluded by the relic density constraint, since the mass difference 
between the two doublet-dominated neutral fermion states is relatively small in this case, leading to efficient 
Xi~X^ annihilation. However, viable parameter points are obtained for larger Yukawa couplings, 1 ^ J/ < 3, 
where the upper bound is imposed to ensure perturbativity, and for masses < 220 GeV. 

As already mentioned in Sec. |lV[ the DM-nucleon cross-section is loop-suppressed in this scenario. Neverthe¬ 
less, due to the large values of the Yukawa coupling y in the viable parameter space, there are still important 
constraints from direct detection experiments. The bounds from LUX |3S] are shown in Fig. With the 
improved sensitivity of XENONIT [35], almost the entire parameter space of the model can be covered. 

At the LHC, the heavier neutral and charged fermions can be produced according to ( [47| , leading to a 
3£ signal with missing energy. Since the mass differences are relatively large in the cosmologically viable 
parameter region, no hard initial-state jet is required. Both ATLAS and CMS have published limits on the 
production cross-section from 3£ searches [151133] (see also additional material available in Ref. [H]), which 
can be compared to the cross-section for the Majorana doublet-triplet model computed with CalcHEP. It is 
found that the current data from ATLAS and CMS does not lead to any limit in the parameter region shown 
in Fig. 

Additionally, this model is constrained by LHC data on the branching fraction of the Higgs boson to two 
photons, as discussed in Sec. jVj Since we require relatively light masses and large Yukawa couplings of the new 
fermions to obtain the correct relic density, the correction to is sizeable. In fact, one finds < 0.6 for the 
whole viable parameter region shown in Fig. (see also Ref. [55]). This is in conflict with results from ATLAS 
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and CMS yielding = 1.17 ± 0.27 [75] and R-y = l-ld^oJa [Sl| the 68% C.L., respectively. Assuming 
that the experimental uncertainties are Gaussian distributed, one thus hnds that the Majorana doublet-triplet 
model with pure doublet dark matter is excluded at more than 95% C.L.. 

However, it is worth pointing out that the h —^ 77 branching fraction may be modified by other new physics 
unrelated to the DM sector, so that the Majorana doublet-triplet model may remain viable as part of a larger 
theory. In this case, the model can be probed more robustly through direct searches for the dark sector 
fermions with future LHC data. For the projected reach of 3£ searches, we use Fig. 10 in Ref. [551 • To account 
for the smaller production cross-section due to the mixing between the dark fermions, we conservatively use 
the 5cr contour in that figure in lieu of the 95% C.L. bound. Using these results, the estimated reach of LHC14 
with 300 fb“^ is illustrated by the short-dashed lines in Fig. For 3000 fb“^, we hnd that the entire viable 
parameter region in the left panel of Fig. [^ can be excluded with 3^ data, while only a small corner in the 
lower left part of the right panel remains. 


VII. CONCLUSIONS 

When probing the fermionic dark matter Higgs portal at colliders, the mediators cannot be integrated out, 
but rather must be considered as dynamical degrees of freedom. The inclusion of new signatures due to the on- 
shell production of the mediators can substantially extend the reach of collider searches for Higgs portal models. 
Using the example of three simple UV completions of the Higgs portal with scalar and fermion mediators, we 
have demonstrated the complementarity of direct detection experiments and high-energy colliders. 

In the case of a scalar mediator, there is no direct connection between direct detection and collider observ¬ 
ables, due to a large number of free parameters. For a singlet scalar, the observed thermal relic density can 
be accommodated with small Higgs-singlet mixing, especially if DM is heavier than the heavy scalar. Small 
mixing helps to evade direct detection bounds, as well as collider bounds from Higgs coupling measurements 
and direct searches for the heavier scalar. Future colliders, i.e. the HL-LHC, ILC or FCC-hh, and direct 
detection experiments can cover a large part of the model’s parameter space. They will, however, not be able 
to test it conclusively. 

In models with fermion mediators, obtaining the correct relic density requires either large Yukawa couplings 
or co-annihilation. The scenario with large Yukawa couplings is generally strongly constrained by direct 
detection, as are models with unsuppressed couplings to the Z boson. The co-annihilation scenario is viable 
for small Yukawa couplings y < 0.1. Since it implies small mass splittings among the lightest states, dedicated 
search strategies are needed at colliders. These compressed spectra can be searched for in mono-jet signatures 
dressed with soft leptons, which are produced in the decay chain of the mediators. 

For singlet DM and a doublet fermion mediator, several scenarios can satisfy all current bounds. If DM is a 
Majorana fermion, Z-mediated DM-nuclei interactions are absent and the co-annihilation scenario is viable for 
DM-mediator mass differences of about 10 — 30 GeV. While the LHC can access DM masses below 300 GeV, 
the scenario can be conclusively tested with XENONIT and future colliders up to TeV-scale DM masses. In 
the Higgs resonance region, rriy^ ~ mh/‘2', the correct relic density can be obtained without co-annihilation for 
Yukawa couplings y ^ O.I, if mediators are below the TeV scale. The LHC (and future hadron colliders) can 
probe mediators up to (beyond) the TeV scale with three-lepton signatures, leaving only a small corner of 
parameter space that neither colliders nor XENONIT can reach. 

If, in turn, DM is a Dirac fermion, the singlet-doublet model is viable only for tiny Yukawa couplings y < 
10“^, where Z-mediated DM-nuclei interactions are sufficiently suppressed to evade LUX bounds. XENONIT 
can probe even smaller couplings. The LHC will test the region < 200 GeV for arbitrarily small Yukawa 
couplings. The ILG and ECC-hh can cover higher masses and ultimately test the model conclusively. 

In a model with one doublet and one triplet fermion held, which both carry weak charges, the only viable 
scenario is pure-doublet DM, where Z- and Higgs-mediated DM-nuclei interactions are absent at tree level. 
Despite the suppression, loop-induced scattering still leads to signihcant constraints from direct detection 
searches. In addition, the model in its minimal version is already excluded by the present Higgs property 
measurements at the LHC, due to the contributions of charged fermions to the Higgs decay to two photons. 
If these can be canceled by a new contribution unrelated to the dark sector, the extended model is viable for 
large Yukawa couplings y > 1 and DM masses below 200 GeV. XENONIT and direct searches for electroweak 
particles at LHC 14 can cover almost the entire viable parameter space of this model. 

Since the mediators in our models are generally not much heavier than the DM particles, an effective Higgs 
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2cw 

-|-i|cos6>2 -|-^sin(202) 

+ ;^sin(26»2) 

0 

+ cos( 26»2) 


TABLE II: Dark fermion couplings to SM bosons in the Majorana doublet-triplet model. 


portal description is not useful. Exceptions are the Higgs resonance region in the Majorana singlet-doublet 
model and parts of the parameter space in the singlet-singlet model, where the DM-mediator mass splitting 
can be sizeable. 

In summary, collider searches and direct detection experiments can pin down models with scalar mediators 
that have moderate mixing with the Higgs boson, and can fully explore models with fermion mediators, apart 
from small regions of parameter space. This is possible due to the combination of different experiments. The 
complementarity of probes will become even more important in the case of a discovery. 
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Appendix A: Loop-induced Higgs boson interactions with Majorana DM 


In this appendix, we give analytic formulae for the loop-in duce d Higgs coupling to DM in the Majorana 
doublet-triplet model with doublet-like DM, Xe = Xb HC). The relevant term in the Lagrangian can 

be written as 


-y AbXb^ + h.c.. 


(Al) 


We have calculated the Yukawa correction Sy at the one-loop level in unitarity gauge in the limit of zero 
momentum transfer to the Higgs boson. The result is given by 


Sy= T. Xj) + 2 gfb^Fiw.xt.xt) 

i,j=a,c 

where the couplings g etc. are listed in Tab. The loop functions can be expressed as 


(A2) 


FiB,Xi,Xj) = 


1 1 \ rrii + m-j 
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p.R R ^ _Lj + m^mp - mf - 2 m| , ^2 n 
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'■B 


mo (mi + mR)2 —m| 


(A4) 
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+ mj{m‘f + 2mimD — ml,) + m%{mYi + 2mimD — m|) — 
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where mu = m° is the doublet DM mass. The notation for the loop integrals Aq and Bq has been adopted 
from Ref. m- Since Xb self-energy corrections do not contribute to the hx^Xt vertex at one-loop, the result 
for Sy as given above is finite without renormalization. 

We have verified that our result is gauge invariant by calculating 6 y in gauge and expanding the result 
in powers of the Higgs mass m^. The leading term Sy{mh = 0) in gauge corresponds to dy in unitarity 
gauge. We have checked that the remainder, 6 y{mh) — Sy{mh = 0), is numerically subleadin^n the parameter 
space where Xb is the lightest neutral state and satisfies the relic density constraint (see Fig. This indicates 
that electroweak box diagrams, which contribute to the effective scalar interaction in Eq. (37) at the same 


order of mb as this remainder, are suppressed with respect to the leading Higgs vertex contributions, which 
enter oc in the amplitude. Our calculation has been performed with two independent computer codes, 
one of which is based on the programs FeynArts m and FeynCalc m- Our result agrees with Ref. [26] . 
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